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Abstract 

For nonnegative real numbers a, (3, 7, A, B and C such that i? + C > and a + /3 + 7 > 0, the 
difference equation 

a + /3x„ + 7a;„_i 

^"+1 /I I p — r7=i ' ■^ = 0,1,2,... 

has a unique positive equihbriuni. A proof is given here for the following statements: 

Theorem 1. For every choice oj positive parameters a, p, 7, A, B and C, all solutions to the difference 
equation 

Xn+i = . , D r-^ > n = 0,l,2,..., a;o £ [0, cxj) 

A + BXn + C Xn-l 

converge to the positive equilibrium or to a prime period-two solution. 

Theorem 2. For every choice oJ positive parameters a, [3, 7, A, B and C, all solutions to the difference 
equation 

a;„+i = — — , 71 = 0,1,2,..., x_i,a;o e (0,oo) 

converge to the positive equilibrium or to a prime period-two solution. 
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1 Introduction and Main Results 



In their book [llj, M. Kulenovic and G. Ladas initiated a systematic study of the difference equation 

a + /3x„ + 7x„_i . , 

^"+1 = 4 I R TT^ ' "- = 0,1,2,... (1) 

A + B Xn + C Xn-l 

for nonnegative real numbers a, (3, 7, A, B and C such that B+C > and a+(3+j > 0, and for nonnegative 
or positive initial conditions x-i, xq. Under these conditions, ([T]) has a unique positive equilibrium. One 
of their main ideas in this undertaking was to make the task more manageable by considering separate 
cases when one or more of the parameters in ([T]) is zero. The need for this strategy is made apparent by 
cases such as the well known Lyness Equation [l7], |22j . [^. 

O -\- Xfi , , 

Xn+l = (2) 

^n— 1 

whose dynamics differ significantly from other equations in this class. There are a total of 42 cases that 
arise from ([T]) in the manner just discussed, under the hypotheses i? + C > and a + /3 + 7 > 0. The 
recent publications [1], [2] give a detailed account of the progress up to 2007 in the study of dynamics of 
the class of equations ([T]). After a sustained effort by many researchers (for extensive references, see [1], 
[2]), there are some cases that have resisted a complete analysis. We list them below in normalized form, 
as presented in |lj , |2j . 

a + Xn 
-\- C Xi-i—i 

a + (3 Xn + J Xn-l 
Xn+1 = (5) 

- a + ^n /„N 

"^'^^^ - A + B Xn + Xn-1 ^ ' 

P Xfi ~\~ Xfi — \ 

"^"^^ - A+BXn + Xn-l 

a + l3Xn + 1Xn-l 
^n+l = -7-— (8) 



A + x 



n—l 



a + Xn+'JXn-l 

Xn+1 = ^ (9) 

O. -\- f3 Xfl H~ X^i—l y V 

^"+1 = ^ , , ^ (10) 



j4 + i? + x< 



n-1 



The dynamics of Equation ([T]) has been settled recently in |3], [20]. Global attractivity of the positive 
equilibrium of Equation ^ has been proved recently in [21]. Since Eq.([6| can be reduced to Eq.Q 
through a change of variables [12], global behavior of solutions to ([S]) is also settled. Equation ([5]) is 
another equation that can be reduced to (|3| , through the change of variables x„ = y„ + 7 [TU] . 

Ladas and co-workers [T], [2], have posed a series of conjectures on these equations. One of them 
is the following. 

Conjecture [Ladas et aL] For equations ^ and (10), every solution converges to the positive equilibrium 
or to a prime period-two solution. 

In this article, we prove this conjecture. Our main results are the following. 
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Theorem 1 For every choice of positive parameters a, (3, 7, A, B and C , all solutions to the difference 
equation 

a + (3 Xn + Xn-l 

= — — — , n = 0,1,2,..., x_i,xo G 0,cx)) (3-3) 

A + B Xn + C Xn-l 

converge to the positive equilibrium or to a prime period-two solution. 



Theorem 2 For every choice of positive parameters a, (3, 7, A, B and C , all solutions to the difference 
equation 

Xn+i = —^ , n = 0,1,2,..., x_i,xo G (0,00) (11) 

D Xn + O Xn—\ 

converge to the positive equilibrium or to a prime period-two solution. 



A reduction of the number of parameters of Eq.( 11 ) is obtained with the change of variables Xn = -^jV'. 



which yields the equation 



yn+i = , , n = 0,1,2,..., y_.i,yo e (0,cx)) (3-2) 

qyn + yn-i 

where r = p = ^, and q = ^. 

The number of parameters of Eq.(3-3) can also be reduced, which we proceed to do next. Consider 
the following affine change of variables which is helpful to reduce number of parameters and simplify 
calculations: 

''-=(^^bTc) y-'BTc^ 



With (12), Eqn.(3-3) may now be rewritten as 



yn+i = , , n = 0,l,2,..., y-i,yo e[L,oo) (3-2-L) 

Qyn + Vn-l 



where 

_ AB+{B+C)I3 
f ~ ylC+(B+C)7 



a - ^ 

q — c 



L 



C(B+C){B a+C a- A 13- A 7) 
(7{B+C)+AC)2 ' 

AC 



(13) 



AC+iB+C)-^ 

Theorems [T] and [2] can be reformulated in terms of the parameters p, q and r as follows. 



Theorem 3 Let a, (3, 7, A, B and C be positive numbers, and letp, q, r and L be given by relations (13). 
Then every solution to Eqn. (3-2-L) converges to the unique equilibrium or to a prime period-two solution. 

Theorem 4 Letp, q, r be positive numbers. Then every solution to Eqn. (3-2) converges to the unique 
equilibrium or to a prime period-two solution. 
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In this paper we prove Theorems |3] and |4j Theorems [T] and |2] follow as an immediate corollary. 

The two main differences between Eq.(3-2-L) and Eq.(3-2) are the set of initial conditions, and the 
possibility of having a negative value of r in Eq.(3-2-L), while only positive values of r are allowed in 
Eq.(3-2). Nevertheless, for both Eq.(3-2-L) and Eq.(3-2) the unique equilibrium has the formula: 

p+l + A/(p+l)2+4r(g + l) 

V = 



2(g + l) 

Although it is not possible to prove Theorem [T] as a simple corollary to Theorem [2] the changes of variables 
leading to Theorems [3] and |4] will result in proofs to the former theorems that are greatly simplified. 

Our main results Theorem [T] and Theorem |2] imply that when prime period- two solutions to Eq.(3-3) 
or Eq.(3-2) do not exist, then the unique equilibrium is a global attractor. We have not treated here 
certain questions about the global dynamics of Eq.(3-3) and Eq.(3-2), such as the character of the prime 
period-two solutions to either equation, or even for more general rational second order equations, when 
such solutions exist. This matter will be treated in an upcoming article of the authors [3]. 

This work is organized as follows. The main results are stated in Section [T| Results from the literature 
which are used here are given in Section [2] for convenience. In Section |3] it is shown that either every 
solution to Eq.(3-2-L) converges to the equilibrium, or there exists an invariant and attracting interval 
/ with the property that the function /(x, y) associated with the difference equation is coordinate-wise 
strictly-monotonic on I x /. In Section |4| a global convergence result is obtained for Eq.(3-2) over a 
specific range of parameters and for initial conditions in an invariant compact interval. Theorem [3] is 
proved in Section |5] and the proof of Theorem |4] is given in Section [6| Section [7] includes computer algebra 
system code for performing certain calculations that involve polynomials with a large number of terms 
(over 365,000 in one case). These computer calculations are used to support certain statements in Section 
|4) Finally, we refer the reader to [IT] for terminology and definitions that concern difference equations. 

2 Results from the literature 

The results in this subsection are from the literature, and they are given here for easy reference. The first 
result is a reformulation of Theorems (1.4.5) — (1.4.8) in |11] . 

Theorem 5 ( [13] . Suppose a continuous function f : [a, 6]^ — > [a,b] satisfies one of i.~iv.: 

i. f{x, y) is nondecreasing in x, y, and 

V(m, M) £ [a, bf, ( /(m, m) = m & /(M, M)= M) m = M 

a. f{x, y) is nonincreasing in x, y, and 

V(m, M) E [a, 6]^ ( /(m, m) = M & /(M, M) = m) =^ m = M 

Hi. f{x, y) is nonincreasing in x and nondecreasing in y, and 

V(m, M) G [a, hf, ( /(m, M) = M k f{M, m) = m) =^ m = M 

iv. f{x, y) is nondecreasing in x and nonincreasing in y, and 

V(m, M) e [a, 6]^ ( /(M, m) = M k f{m, M) = m) =^ m = M 
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Then Un+i = f{yn,yn-i) has a unique equilibrium in [a,b], and every solution with initial values in [a,b] 
converges to the equilibrium. 

The following result is Theorem A. 0.8 in 

Theorem 6 Suppose a continuous function f : [a, 6]^ [a, 6] is nonincreasing in all variables, and 

V(m,M) G [a,6]^ ( /(m, m, m) = M &; f{m, m, m) = M ) =^ m = M 

Then Un+i = f{yn,yn-i,yn-2) has a unique equilibrium in [a,b], and every solution with initial values 
in [a, b] converges to the equilibrium. 

Theorem 7 ( [5j ) Let I be a set of real numbers and let F : I x I ^ I be a function F{u, v) which decreases 
in u and increases in v. Then for every solution of the equation 

Xn+l = F{Xn,Xn^l), n = 0, 1, . . . (14) 

the subsequences {x2n} o,nd {x2n+i} of even and odd terms do exactly one of the following: 

(i) They are both monotonically increasing. 

(ii) They are both monotonically decreasing. 

(iii) Eventually, one of them is monotonically increasing and the other is monotonically decreasing. 
Theorem [7] has this corollary. 

Corollary 1 (jS]) If I is a compact interval, then every solution of Eq.(14) converges to an equilibrium 
or to a prime period-two solution. 

Theorem 8 ([8j) Assume the following conditions hold: 

(i) /i G C[(0,cx)) X (O,oo),(0,cx))]. 

(ii) h{x, y) is decreasing in x and strictly decreasing in y. 

(iii) X h(x, x) is strictly increasing in x. 

(iv) The equation 

Xn+l = Xnh{Xn,Xn-l), n = 0, 1, . . . (15) 

has a unique positive equilibrium x. 



Then x is a global attractor of all positive solutions of Eq.(15) 



3 Existence of an Invariant And Attracting Interval 

In this section we prove a proposition which is key for later developments. We will need the function 

fi ^ r + px + y 

f{x,y) := ■ , x,ye[L,oo) (16) 

qx + y 

associated to Eq.(3-2-L). 

Proposition 1 At least one of the following statements is true: 



(A) Every solution to (3-2-L) converges to the equilibrium. 

(B) There exist m* , M* with L < m* < M* s.t. 

(i) [m*,M*] is an invariant interval for Eq. (3-2-L), i.e., f{[m*,M*] x [m*,M*]) C [m*,M*]. 

(ii) Every solution to Eq. (3-2-L) eventually enters [m*,M*]. 

(iii) f{x,y) is coordinate-wise strictly monotonic on [m*,M*]^. 

The next lemma states that the function /(•, •) associated to Eq. (3-2-L) is bounded. 
Lemma 1 There exist positive constants C and U such that L < C and 

C<f{x,y)<U, X, ye[L,oo) (17) 



In particular. 
Proof. The function 



f {[C,U] i^M) C [C,U] (18) 



associated to Eq.(3-3) is bounded: 



min{a,/3,7} a + f3x + 'yy max{a,/3,7} 



max{A,B,C} - A + Bx + Cy - m.m{A,B,C]' 



Set £ := ^^"^ ^ ui^IabI}} • '^^^ affine change of coordinates (jl2j) maps the rectangular 

region \L,W^ onto a rectangular region 

{LMf which satisfies (rm and (M. □ 



Lemma 2 If p = q, then every solution to Eq. (3-2-L) converges to the unique equilibrium. 

Proof, p = q then Dif{x,y) = — J_^^^2 and D2f{x,y) = — ■ Thus, depending on the sign of r, 

the function f(x, y) is either nondecreasing in both coordinates, or nonincreasing in both coordinates on 
[L,oo). By Lemma [T| all solutions {yn}'^=-i satisfy yn G for n > 1. A direct algebraic calculation 

may be used to show that all solutions (m, M) S [C,U] of either one of the systems of equations 

M = f{M,M) ( M = f{m,m) 

m = f{m,m) \ m = f{M,M) 

necessarily satisfy m = M. In either case, the hypotheses (i) or (ii) of Theorem [5] are satisfied, and the 
conclusion of the lemma follows. □ 

We will need the following elementary result, which is given here without proof. 
Lemma 3 Suppose q ^ p. The function f{x,y) has continuous partial derivatives on {L,oo)'^, and 

i. Dif{x,y) = if and only if y = and Dif{x,y) > if and only if {p-q)y> qr. 

ii. D2f{x, y) = if and only if x = and D2f{x, y) > if and only if {q — p) x > r. 

We will need to refer to the values Ki and K2 where the partial derivatives of /(x, y) change sign. 
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Definition 1 If p ^ q, set 



K 



q r 



1 • = 



and K2 



p-q 



Definition 2 For L <m< M, let 



(f){m,M) := mm{f{x,y) : {x,y) G [m,M]^} and ^{m,M) := max{f{x,y) : (x,y) G [m, M]^} 

Lemma 4 Suppose p ^ q. If [m, M] C [C,U] is an invariant interval for Eq.(3-2-L) with m < Ki < M 
or m < K2 < M, then m < (f){m, M) or <I>(m, M) < M or m = M = y. 



Proof. By definition of (j) and <I>, m < (p{m,M) and <I>(m, M) < M. Suppose 

m = (t>{m,M) and ^{m,M)=M 



(19) 



Tlie proof will be complete when it is shown that m = M. There are a total of four cases to consider: (a) 
r > and p > q, (h) r < and p < q, (c) r > and p < q, and (d) r < and p > q. We present the proof 
of case (a) only, as the proof of the other cases is similar. 

If r > and p > q, then Ki G [m, M] and K2 [m, M]. Note that 

[m, M] X [m, M] = [m,M] X [m,Ki] |J [m,M] x [/s:i,M]. 

By Lemma [3] the signs of the partial derivatives of f{x, y) are constant on the interior of each of the sets 
[m, M] X [m, i^Ti] and [m, M] x [iTijM], as shown in the diagram. 

(m,M) 



(m, m) 



/(T,i) 

/a,i) 



Figure 1: The arrows indicate type of coordinate- wise monotonicity of /(x, y) on each region. 
Since /(x, y) is nonincreasing in both x and y on [m, M] x [m, ETi], 

f{M,Ki) < f{x,y) < f{m,m) for G [m, M] x 

Similarly, f(x,y) is non decreasing in x and nonincreasing in y on [m, M] x [Xi,M], hence 
f{m,M) < f{x,y) < f{M,Ki) for {x,y) G [m, M] x [Ki,M]. 



From (20) and (21) one has 



{m,M)=f{m,M) and ^{m, M) = f{m,m) 



Combine (22) with relation (19) to obtain the system of equations 

/(m, M) = m 

f{m,m) = M 



(20) 



(21) 



(22) 



(23) 
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Eliminating M from system (23) gives the cubic in m 

g (g + 1) + (1 — J"?) fn^ + (—1 — p — qr)m — r = (24) 



which has the roots 



1 l-p-\/(l+p)^ + 4r(l + g) 1 -P+ +4r(l + 

, 3,110. 



2(1 + 



2(l + g) 



(25) 



Only one root in the list (25) is positive, namely 

1-P + V(l +p)2 +4r(l + 



m 



2(1 + 



Substituting into one of the equations of system ( |23[ ) one also obtains M = y, which gives the desired 
relation m = M = y. □ 



Definition 3 Let mo := C, Mq :=l(, and for £ = 0,1,2,... let rrii := (j){mi,Mi), := ^>(m<?,M£). 

By the definitions of m£. Mi, •) and <!>(•, •), we have that [m^+i, M^j^i] C [m^, M^] for £ = 0, 1, 2, . . .. 
Thus the sequence {mi} is nondecreasing and {M{\ is nonincreasing. Let m* := limm^ and M* := limM^. 



Lemma 5 Suppose p ^ q. Either there exists N £N such that {Ki, K2}ri[mN, M^] 



or m 



M* = y. 



Proof. Arguing by contradiction, suppose m* < M* and for all £ G N, {Ki,K2} n [mi,M£\ ^ 0. Since 
the intervals \mi,M^ are nested and ^\mi,M^ = [m*,M*], it follows that {^1,^2} n [m*,M*] / 0. By 
Lemma |4j we have 

m*<4>{m*,M*) or $(m*,M*)<M* (26) 
Continuity of the functions </> and ^ implies 

(p{m*,M*) = lim (l){mi. Mi) = limm^+i = m* or ^>(m*,M*) = lim <I>(m^, M<?) = limA%i = M* (27) 



Statements (26) and (27) give a contradiction. 



□ 



Proof of Proposition [TJ Suppose statement (A) is not true. By Lemma [2j one must have p ^ q. Note 
that if {yi} is a solution to Eq.(3-2-L), then yg^i G [mi. Mi] for £ = 0, 1, 2, . . .. If m* = M* , since mi — > m* 
and Ml M* we have yi — > y, but this is statement (A) which we are negating. Thus m* < M* , and by 
Lemma [4] there exists G N such that {-f^i, K2} H [m^r, M]\f] = 0, so f{x, y) is coordinate- wise monotonic 
on [m7v,M7v]. The set [mp^,Mp^] is invariant, and every solution enters [mN,Mi\j] starting at least with 
the term with subindex A + 1. We have shown that if statement (A) is not true, then statement (B) is 
necessarily true. This completes the proof of the proposition. □ 
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4 The Equation (3-2) with r>0,p>q and ^ < f 

In this section we restrict our attention to the equation 

Vn+i = f{xn,Xn-i) n = 0, 1, . . . , x_i , G (0, oo) (3-2) 

where 

f( N r+px + y 

j{x,y) := , 

qx + y 

For p > 0, q > 0, and r > 0, Eq.(3-2) has a unique positive equiUbrium 

_ P+1+ y^{p+l)^+4r{q + l) 

' = Hi[TY) 

We note that if I C (0, oo) is an invariant compact interval, then necessarily y €z I. 

The goal in this section is to prove the following proposition, which will provide an important part of 
the proofs of Theorems [T] and |4j 

Proposition 2 Let p, q and r be real numbers such that 

p>q>0, r > 0, and < - (29) 

p-q q 

and let [rh,M] C (^^)|) be a compact invariant interval for Eq. (3-2). Then every solution to Eq. (3-2) 
with X-i,xq G [m, M] converges to the equilibrium. 

Proposition |2] follows from Lemmas |6j [7] and [9j which are stated and proved next. 



Lemma 6 Assume the hypotheses to Proposition [1[ If either q > 1 or p < 1, then every solution to 
Eq. (3-2) with X-i,xq £ [rh,M] converges to the equilibrium. 

Proof. We verify that hypothesis (iv) of Theorem 5 is true. Since x > for x G [m, M], the function 



(30) 



f(x,y) is increasing in x and decreasing in y for (x,y) G [m,M]^ by Lemma [sj Let m, M G [m, M] be 
such that m ^ M and 

f(M,m) - M = 

f(m,M)-m = 



We show first that system (30) has no solutions if either q > 1 oi p < 1. By eliminating denominators in 

° (31) 



both equations in (30), 



m-mM + Mp - M'^q + r 
M — mM + mp — m? q + r 



and by subtracting terms in (31) one obtains 

(M - m)(l - p + g (m + M)) 







(32) 



Since m ^ M, we have q (m + M) = p — 1, which implies that for p < 1 there are no solutions to system 
(30) which have both coordinates positive. Now assume p > 1; from (32), m = En^^d^Ll^ and substitute 
tne latter into ( 30 ) to see that x = M is a solution to the quadratic equation 



(1-q) x' + tl±Aizl±A^ + zl±p±ll = o 



(33) 
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By a symmetry argument, one has that x = m is also a solution to (33). By inspection of the coefficients 



of the polynomial in the left-hand-side of ( 33 ) one sees that two positive solutions are possible only when 



q < 1. To get the conclusion of the lemma, note that the fact that (30) has no solutions with m 7^ M is 



just hypothesis (iv) of Theorem [5j □ 

Lemma 7 Assume the hypotheses to Proposition^ If 

r < q — p (34) 
then every solution to Eq. (3-2) with X-i,xq G [m, M] converges to the equilibrium. 
Proof. By substituting Xn = /(x„_i,x„_2) into Xn+i = /(x„,x„_i) we obtain 

I r+pXn-l+Xn-2 I 

_ r+pXn + Xn-l _ qx„_i+x„_2 + ^"-1 _ 

^ri+l- , „ ^ - r+px„_i+x„_2 ' n-U,i,... {6ii) 

qXn + Xn-l q +Xn-l 

that is, 

?/ . . pr+p^y + qry + qy'^ +pz + rz + yz 

Xn+i = f{xn,Xn-i,Xn-2), where / (x, y, z) = ■ ■ 5- ■ (36) 

qr+pqy + qx'^ + qz + yz 

where the x has been kept in /(x, y, z) for bookkeeping purposes. Thus /(x, y, z) is constant in x. We 
claim /(x, y, z) is decreasing in both y and z. To see that the partial derivative 

D-J{x,y,z) = Jl±(l^imzll±Sl^ (37) 
[qr + pqy + qy"^ + q z + y z) 

is negative just use p > q and the inequality {p — q) y — qr > 0, which is true by Lemma [sj The remaining 
partial derivative is 

L{y,z) 

{qr + pqy + qy"^ + qz + y z) 

where 



D2KX, y, z) = ■ ■ ^ (38) 



/i(y, z) := —q^r'^ + 2pqry — 2q^ry+p'^qy'^—pq'^y'^ + q^ry'^+prz — qrz + 

We have. 



ry ry O 'J O 

+ pqrz — qrz + 2pqyz — 2qyz + 2qryz+pz— qz+rz 



Dih{y, z) = 2{p — q)qr + 2q {p^ — pq + qr) y + 2q{p — q + r) z > 0, (39) 
D2h{y,z) = {p-q) {1 + q) r + 2q {p-q + r) y + 2 {p-q + r) z>0 (40) 



ir qr \ q{l + qf r"^ {p^ - p q + qr) 



Since < m, 

p-q — ' 

h{y, z)>h[ \ = " ' ^ • ^ > , y, zG[m,M] 

Vp-g p-qj {p-q) 

thus we conclude that D2f{x, y,z) < for x,y, z £ [m, M]. 

To complete the proof we verify the hypotheses of Theorem |6j We claim that the system of equations 

f{m,m,m)-M = , . 

f{M, M,M)-m = ^ ' 
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has no solutions (m,M) with m ^ M whenever hypothesis (34) holds. 
By eliminating denominators in both equations in (41) one obtains 

--m^ + M — mp — mp'^ — q + m AI q + AI q + m M pq — mr — pr — rnqr + M qr = 
—AP + m A'P — M p — AIp^ + mAI q — AP q + mAP q + m M p q — Al r — pr + mqr — M qr = 



and by subtracting terms in ( 42 ) one obtains 

((m — M) (— m — M + m M — p — — m q — M q + m M q — r — 2 qr)^ = 



(42) 



(43) 



Since m ^ M, we may use the second factor in the left-hand-side term of (43) to solve for M in terms of 
m, which upon substitution into /(m, m, m) = M and simplification yields the equation 



a2 m + aim + uq 



( — 1+m) (1 + q) {m? + mq + m? q + mpq + qr) 







(44) 



where 



ai 

02 



p + 2pq + p q + qr + 2q rj 



p + p + 2pq + 3p q + p q + r — pr + 4qr + 4q r + 2pq r 



{1 + q) {l + 2q + pq + qr) 



By hypothesis (34) we have rp<pq—p < p q, hence p q — rp > 0, which implies ai > 0. By direct 
inspection one can see that > and 02 > 0. Thus (44) has no positive solutions, and we conclude that 
(41) has no solutions (m, M) G [m, M] with m 7^ M. We have verified the hypotheses of Theorem [sj and 
the conclusion of the lemma follows. □ 



Lemma 8 Let p > 0, q > and r >0. If the positive equilibrium y of Eq.(3-2) satisfies y < ^, then y is 
locally asymptotically stable (L.A.S.). 

Proof. Solving for r in 

_ _ r + {p+l)y 

gives 

r = {q+l)y^-{p + l)y (45) 

Then a calculation shows 

Dif{y,y) = ^ 
D2f{y,y) = -yf^) 

Set ti := Dif{y,y) and t2 := D2f{y,y). The equilibrium y is locally asymptotically stable if the roots of 
the characteristic polynomial 

p{x) = — tiX — t2 

have modulus less than one [11]. By the Schur-Cohn Theorem, y is L.A.S. if and only if < 1 — ^2 < 2. 
It can be easily verified that 1 — i2 < 2 if and only if < q y + 1 which is true regardless of the allowable 
parameter values. Since p — > by the hypothesis, we have = | ^{^qli) I ~ ^^+1) ' ^^^^^^ some algebra 
gives < 1 — t2 if and only if 

Ip+l 



2, + l<5 

But (46) is a true statement by formula (28). We conclude y is L.A.S. □ 
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Lemma 9 Assume the hypotheses to Proposition^ If 

p > 1, 9 < 1 , and r > p^ q — p 
then every solution to Eq. {3-2) with X-i,xo £ [m,M] converges to the equilibrium. 



(47) 



Proof. The proof begins with a change of variable in Eq.(3-2) to produce a transformed equation with 
normalized coefficients analogous to those in the standard normalized Lyness' Equation [T7], [22], [9] 



^n—l 



(48) 



We seek to use an argument of proof similar to the one used in pT] , in which one takes advantage of the 
existence of invariant curves of Lyness' Equation to produce a Lyapunov-like function for Eq.(3-2). 
Set yn = pZfi in Eq.(3-2) to obain the equation 



a + Zn + gZn-1 



bZn + Z, 



n—1 



where 



p 



n = 0,1, 
1 

9= -, 
p 



z-i, zo G (0, oo) 



We shall denote with z the unique equilibrium of Eq.(49). Note that 



z = py 



(49) 
(50) 

(51) 



It is convenient to parametrize Eq.(49) in terms of the equilibrium. We will use the symbol u to represent 
the equilibrium z of Eq.(49). By direct substitution of the equilibrium u = z into Eq.(49) we obtain 

r = {b + l)u'^ - {g + l)u (52) 



By (|52[, r > iff ti > f±i. Using (|52 

9+1 
6+1 

{b+l)u'^ - {g + l)u + Zn+g Zn-l 



g > and u > equivalent to Eq.(|49|): 



to eliminate r from Eq.(49) gives the following equation for 6 > 0, 



b Zfi -\- Zyi—\ 



n = 0,1,2, 



yo G (0,oo) 



(53) 



Therefore it suffices to prove that all solutions of Eq.(53) converge to the equilibrium u. 
The following statement is crucial for the proof of the proposition. 

Claim 1 u > 1 if and only if r > p^ q — p. 

Proof. Since y = p^ = pu, we have u > 1 if and only if y > which holds if and only if 

P+1 + \/(p+l)2+4r(g + l) 

WTT, >" 

After an elementary simplification, the latter inequality can be rewritten as r > p^ q — p. 
By the hypotheses of the lemma, by Claim [T] and by (50) and (52) we have 



□ 



6<1, g <l, 1 <u < 



and < u 

b+1 - 



(54) 
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We now introduce a function which is the invariant function for (48) with constant a. 



case the the equihbrium of (48) is u): 



(l + x)(l+y)(n2-M + j; + ?/) 



xy 



u (in this 



(55) 



Note that g{x,y) > for all x,y £ (0,oo) whenever u > 1. By using elementary calculus, one can show 
that the function g{x,y) has a strict global minimum at {u,u) |9], |22], i.e., 



g{u,u) < g{x,y), (x, y) G (0, oo)^ 
We need some elementary properties of the sublevel sets 

S{c) := {(s, t) E (0, oo) : g{s, t) < c} , c> 
We denote with Qi{u,u), i = 1,2,3,4 the four regions 



Qi{u,u) 
Q2{u,u) 
Qz{u,u) 
Qa{u,u) 



{{x,y) G (0,oo) X (0,oo) 
{(x,y) G (0,oo) X (0,oo) 
{{x,y) G (0,cx)) X (0,cx)) 
{{x,y) £ (0,cx)) X (0,oo) 



Let 



T{x,y) := { y, 



{g+l)u^ - {b+l)u + y + gx 



by + X 

be the map associated to Eq. ( [SS] ) (see [16]). 

Claim 2 If {x,y) £ Q2{u,u) U Qi{u,u) \{{u,u)}, then g{T{x,y)) < g{x,y). 
Proof. Set 

Ai(x, y) := g{x, y) - g{T{x, y)) 



u < X, u < y } 

x < u, u < y } 

X < u, y < u } 

u < X, y < u } 



{x,y) G (0,oo) X (0,oo) 



A calculation yields 



Ai(x,y) 



jl + x) Fi{x,y)F2{x,y) 
xy {hx + y) F3{x,y) 



where 



Fi{x,y) 
F2{x,y) 
F3{x,y) 



h{x - u) {y - \) + {y - u) {hu + y + u - g) 

b{x — u)'^ + b {x — u) u + b {x — u) + {u — y){b + y g) 

{b+l)u^ -{l+g)u + x + gy 



(56) 



(57) 



(58) 
(59) 



By (54), for {x,y) G Qi{u,u) \ {{u,u)} we have u < x and y < u < ^ with {x,y) ^ {u,u), therefore 



Fi{x~y) < 0, F2{x,y) > and F^^x^y) > 0. Consequently Ai{x,y) > for {x,y) G Qi{u,u) \ {(u, n)}. To 
see that Ai(x,y) > for {x,y) G Q2iu,u) \ {{u,u)} as well, rewrite Fi{x,y) and F2{x,y) as follows: 

Fi{x,y) = b{u-x){\-u) + {y-uf + {y-u)u + {y-u){u-g) + b{y-u)x 
F2{x,y) = b{x - u){x + u^) -b{y - u)u^ - {y - uf g - {y - u)ug 



For (x, y) G Q2{u, u) \ {(n, u)} we have x < u < y and (x, y) / (n, u). Thus Fi{x, y) > 0, F2{x, y) < 0, and 
F?,{x,y) > 0, which imply Ai(x,y) > 0. □ 
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Claim 3 Suppose g > b. If (x, y) G Qi{u, u) U Qz{u, u) \ {{u, u)}, then g{T'^{x, y)) < g{x, y). 

Proof. This proof requires extensive use of a computer algebra system to verify certain inequalities 
involving rational expressions. Here we give an outline of the steps, and refer the reader to Section [7] for 
the details. 

Since b<g<\<u<j^, and < u we may write 

u = +t, i > 

(60) 

9 = Si' 

The expression A2 := g{x,y) — g{T^{x,y)) may be written as a single ratio of polynomials, A2 = ^ with 
D > 0. The next step is to show > for (x, y) 7^ {u, u). 

Points (x,y) in Qi{u,u) may be written in the form x = u + v, y = u + w, where v,w £ [0,oo). 
Substituting x, y, u and g in terms of v, w, s and t into the expression for one obtains a rational 
expression ^ with positive denominator. The numerator N has some negative coefficients. At this points 
two cases are considered, w > v, and w < v. These can be written as w = v + k and v = w + k for 
nonnegative k. Substitution of each one of the latter expressions in gives a polynomial with positive 
coefficients. This proves A2(x,y) > for (x,y) S Q2{u,u). 

If now we assume (x, y) G Qs^u, u) with (x, y) 7^ (u, u), we may write 

X = ^ t;G[0,oo) 

(61) 

y = ^ w£[0,oo) 

The rest of the proof is as in the first case already discussed. Details can be found in Section [7} □ 



Claim 4 Suppose g < b and u > 1. If {x,y) £ Qi{u, u) U Qz{u, u) \ {(u, u)}, then g{T^{x, y)) < g{x, y). 

Proof. The proof is analogous to the proof of Claim |3j We provide an outline. More details can be found 
in Section [7l 

Since 1 < ti, we may write u = 1 + t with f > 0. Also, u < \ implies b < ^, and b = jip^ for s > 0. 
Since g < b we may write g = for £ > 0. 

The expression A3 := g{x,y) — g{T^{x,y)) may be written as a single ratio of polynomials, A2 = ^ 
with D > 0. The next step is to show > for {x,y) 7^ {u,u). This is done in a way similar to the 
procedure described in in Claim [3} □ 

To complete the proof of the lemma, let {(p,ij) G (0, 00) x (0,oo). Let {yn}n>-i be the solution to 
(53) with initial condition (y_i,yo) = [4>^'^)i a-iid let {r"((/), ■(/;)}n>o be the corresponding orbit of T. The 
following argument is essentially the same as the one found in [211 ; we provided here for convenience. 
Define 

c := liminf 5(r"(0, V)) (62) 

n 

Note that c < 00, which can be shown by applying Claims |2] [3] and |4] repeatedly as needed to obtain a 
nonincreasing subsequence of {g(T"((/), V'))}ri>o that is bounded below by g{u,u). Let {g(T"*((/), ^))}fc>o 
be a subsequence convergent to c. Therefore there exists c > such that 

c/(r"'=((/),V)) < c forall/c>0. 
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that is, 

r"H</',^)) G S{c) := {{s,t) : g{s,t) < c} for A; > 
The set S{c) is closed by continuity of g{x,y). Boundedness of S{c) follows from 

n / / {1 + x){l + y){u^ - u + X + y) f \ ^ cr \ 

0<x,y< =g(^x,y) = c, for {x,y) £ S[c) 

xy 

Thus S{c) is compact, and there exists a convergent subsequence {T"''=« ((/>, ^))}f with limit {x,y). Note 
that 

c= lim g{T^''e{(l),i^))=g{x,y) (63) 
£^00 

We claim that {x,y) = {u,u). If not, then by claims [2][3] and |4| 

mm{g{Tix,y)),g{T\x,y)),g{T\x,y))} < c (64) 
Let II • II denotes the euclidean norm. By ( |64[ ) and continuity, there exists d > such that 

\\{s,t)-{x,y)\\<^ =^ mm{g{T{s,t)),g{T\s,t)),g{T'{s,t))} < c (65) 
Choose L E N large enough so that 

\\T'"^L{cp,^)-{x,y)\\<S (66) 



But then (65) and (66) imply 

mm{g{T^>'L+\cf>,i;)),g{T^''L+\s,t)),g{r^>=L+-\s,t))} < c (67) 



which contradicts the definition (62) of c. We conclude {x,y) = {u,u). From this and the definition of 
convergence of sequences we have that for every e > there exists L G N such that \\T"'''l ((f>,'tp) — {u,u)\\ < e. 
Finally, since 



max 



\\ynk^-i - u\,\yn^^ -n|| < ||(ynfe^_i -n,y„^^ -'u)|| = ||r"''-i -0) - {u,u) 



we have that for every e > there exists L G N such that lyn^^ — u\ < e and \ynf,^-i — u\ < e. Since n is a 
locally asymptotically stable equilibrium for Eq.(53) by Lemma [s] it follows that 2/„ u. This completes 
the proof of the lemma. □ 



5 Proof of Theorem [3] 

To prove Theorem [s] it is enough to assume statement (B) of Proposition [T| Also by Lemma [2] we may 
assume p ^ q without loss of generality. Thus we make the following standing assumption, valid throughout 
the rest of this section for Eq.(3-2-L). 

Standing Assumption (SA) Assume p ^ q and that there exist m* , M* with C < m* < M* < lA such 
that for Eq.(3-2-L) and its associated function f{x,y), 

(i) [m*,M*] is an invariant interval. 

(ii) every solution eventually enters [m*,M*]. 

(iii) f{x,y) is coordinate-wise strictly monotonic on [m*,M*]^. 
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The function f{x,y) is assumed to be coordinate- wise monotonic on [m*,M*], and there are four 
possible cases in which this can happen: (a) f(x, y) is increasing in both variables, (b) f(x, y) is decreasing 
in both variables, (c) /(x, y) is decreasing in x and increasing in y, and (d) /(x, y) is increasing in x and 
decreasing in y. 

We present several lemmas before completing the proof of Theorem |3j 
By considering the restriction of the map T of Eq.(3-2-L) to [m*, Af*]^, an application of the Schauder 
Fixed Point Theorem [7j gives that [m*,M*]^ contains the fixed point of T, namely {y,y). Thus we have 
the following result. 

Lemma 10 y e [m*,M*]. 

Lemma 11 Neither one of the systems of equations 

( M = f{M,M) ( M = f{m,m) 

^^^> \ m = f{m,m) """^ ^^^> \ m = /(M,M) 

have solutions {m,M) G [m*,M*]^ with m < M. 

Proof. Since x = y is the only solution to f{x,x) = x, it is clear that only {y,y) satisfies (Si). Now let 
(m, M) be a solution to (S2). From straightforward algebra applied to M — m = /(m, m) — f{M, M) one 
arrives at (p + 1)(M — m) = 0, which implies m = M. □ 

Lemma 12 Suppose f{x, y) is increasing in x and decreasing in y for (x, y) € [m*, M*] . Then p — q > 0. 
Proof. By the standing assumption (SA), p ^ q. By Lemma [sj the coordinate- wise monotonicity hypoth- 



esis, and the fact y G [m*, M*\ from Lemma 10 we have 



{p — q)y > qr and {q — p)y < r (68) 



The inequalities in (68) cannot hold simultaneously unless p — q > 0. □ 



Lemma 13 // /(x, y) is increasing in x and decreasing in y for (x,y) G [m*,M*], then f{[m*,M*]^) C 
(1>|)- 

Proof. For (x,y) G [m*,M*], the function / is well defined and is componentwise strictly monotonic on 
the set [x,oo) x [y, cxd). Then, 

f{x,y) < lim f{s,y) = lim + ^ = P 

s^oo q s -\- y q 

(69) 

SI \ T ^\ V r+px + t 

f{x,y) > hm f{x,t) = lim — — = 1 

t~>oo t-*oo qx + t 

□ 

Lemma 14 Let p > 0, q > and r > 0. If f{x,y) is increasing in x and decreasing in y on [m*,M*], 
then 

^ < P (70) 

p — q q 
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Proof. Since y £ [m*,M*] by Lemma 10, we have Di(jj,y) > and D2{y,y) < 0. By Lemma 12 p > q, 
and by Lemma |3] 

{p — q)y>qr and {q — p)y<r (71) 

Then, 

(72) 



y > 



p-q 



In addition, by Lemma [3j 



y = f{y,y) < lim f{s,y) = lim 



r +ps + y p 



s^oo qs + y 



(73) 
□ 



Lemma 15 Suppose f{x,y) is increasing in x and decreasing in y for {x,y) S [m*,M*]. If r < 0, then 
p — q + r > 0. 



Proof. Since Dif{x,y) > for {x,y) G [m*,M*], and by Lemma |3j Lemma 10 and by Lemma 12 



we 



have y > that is. 



v/(p + l)^ + 4r(g+l) > ^^(^+^^ -p-l 



p-q 



(74) 



If the right-hand-side of inequahty (^74| is nonnegative, then, after squaring both sides of ( 74 1 we have 



{p+lf + 4r{q + l) > 



2r{q + l)Y 4r{q+l){p+l) 



p-q 



+ 



p-q 



Further simphfication of ( 75 ) and the hypothesis r < yield 



^ ^ r{q + l) ^p+1 



(75) 



(76) 



{p-q)"^ p-q 

which, after some elementary algebra, implies p — q + r > 0. Now assume the right-hand-side of inequality 



(74) is negative, relation that we may rewrite as 



— r 1 p -|- 1 

< - 

p-q 2 q+1 



(77) 



If I 2±1 < X, then —^ < 1, which gives the conclusion p — q + r > 0. If| 2±1 > x, that is, p > 2 g -|- 1, then 



p — q + r>q + r+ l 



Therefore ifg-|-r-|-l >0 the conclusion of the lemma follows from this and from ( 78 ) . Assume 

q+r+1 < 



From relations (13) we have 

{b + c) (a'^ + b a + a - a P + 2 ab c-r + a c^ + b'^ + b cy"^ + -f'^) 



+ r + l 



c{ac + bj + c^y 



hence assumption (79) and relation (80) imply 

R-- a'^ c'^ ^bc^ a + a - ac^ P + 2abc-f + ac^ + b^j'^ + 2bcj'^ + -f^ <0 



(78) 
(79) 

(80) 

(81) 
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Further algebra gives 

/ 2 , /? I A 2 , bca-f c^a-f , , 2 , 2 , 267^ , ^^7^ , '^^^ ^ n /qo^ 

— it— (— c a + ac7 — CP7 + 07 )=c aH 1 h 7 + c 7 H 1 1 > (82) 

ac a a a ac a 



Since i? < by (81 ), from inequality (|82| we have 

— c^a + ac7 — c/37 + 67^ <0 



Finally, from (13) we have 



p — q + r 



(b + c) (— Q + ac7 — c/37 + 67^) 



c (a c + 67 + 07)^ 



Combining (83) with (84) we obtain p — q + r > 0. 



(83) 



(84) 



□ 



Lemma 16 // r < and f(x,y) is increasing in x and decreasing in y for {x,y) G [m*,M*], then every 
solution converges to the equilibrium. 



Proof. Since p — q + r > Oby Lemma \15\ we have K2 = — < 1 , which together with Lemma [T3| implies 



that [1, 2] is an invariant, attracting compact interval such that f(x,y) is increasing in x and decreasing 
in y on [1, |]^- Since /([I, |]^) C (1,|), we see that every solution to Eq.(3-2-L) eventually enters the 
invariant interval (1, |)- The change of variables 



Vn 



1 + 2 z 

1 + Zn 



or z. 



n — p 



transforms the equation 



r +pyn + yn+i „ . 

yn+i = , , n = 0,l,., 

q y-n + yn+1 



y-i,yo G 1, 



p 



into the equivalent equation 

Zn+l = g{Zn,Zn-l), 71 = 0,1,.. 

where 

g{w,v) : 



2:_i,Zo G (0,00) 
q{l + v){- {q {p- q + r)) + (-p^ +pq- qr) w) 



(85) 



(86) 



(87) 



(1 + w){—q {p — q — qr) + {—p^ + pq + r) v) 

We claim that for w,v ^ (0,oo), (a) g{w,w) is increasing in w, (b) g{w,v)/w is decreasing in w, and (c) 
g{w, v)/w is decreasing in v. Indeed, since p>q, r<Q,p — q + r>Q, and < 2 we have 

/ / -r q^ (I + q) {-p + qf 

— [g[w,w)) = 2 > 

d-u; (— p q -\- -\- r — p^ w + p qw + q^ r w) 



d ( g{w,v) 
dw 



q {-p + q) iq {p - q + r) + ip^ - p q + qr) w) 



< 



d f g{w,v) 
dw \ w 



{—p q + + q^ r — p^ V + p qv + q^ r v)^ w {1 + w) 

q {1 + v) [q {p — q + r) + 2 q {p — q + r) w + {p^ — pq + qr^ u;^) 
{q {p — q — qr) + (p^ — pq — qi f) v) w'^ {1 + w)^ 



< 
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Also, note that Eq.(87) has a unique equilibrium z. Therefore hypotheses (l)-(4) of Theoremplare satisfied, 



so every solution {^n} to Eq.(87) converges to z. By reversing the change of variables, one can conclude 



that every solution to Eq.(86) converges to the equilibrium. 



□ 



Proof of Theorem [3[ The four parts of the proof are: 



a. f{x,y) is increasing in both x and y on [m*,M*]^; By Lemma 11 the hypotheses of Theorem [5] part 
(i). is satisfied, hence every solution converges to the equilibrium y. 



b. /(x, y) is decreasing in both x and y on [m* , M*\ : By Lemma 11 the hypotheses of Theorem^part 



(ii). is satisfied, hence every solution converges to the equilibrium y. 

c. f{x,y) is decreasing in x and increasing in y on [m*,M*]'^: By the corollary to Theorem [t] we 
conclude every solution converges to the unique equilibrium or to a prime period-two solution. 

d. f{x, y) is increasing in x and decreasing in y on [m*,M*]'^: By Lemmas [3 , [l2 and 13 there is no loss 
of generality in assuming [m* , M*] C {K, where K := max{^^, ^^}, which we do. We consider 
two subcases. If r > 0, then Lemma [12] Lemma 14 and Proposition [2] imply that every solution 



converges to the unique equilibrium. If r < 0, then Lemma 16 implies that every solution converges 
to the unique equilibrium. 

This completes the proof of Theorem [3| Since Theorem |3] is just a version of Theorem [T] obtained by an 
affine change of coordinates, we have also proved Theorem [T] as well. □ 

6 Proof of Theorem [4] 

The first lemma guarantees solutions to Eq.(3-2) to be bounded. 

Lemma 17 Let p > 0, q > and r > 0. There exist positive constants C and U such that every solution 
{xn}^_i to Eq.(3-2) satisfies Xn G [>C,Z//] for n>2, and the function 



satisfies 



Proof. Set 



f( \ r+px + y 2 

f{x, y) = , (x, y) G (0, oo) 

qx + y 

f{[CM] X [CM]) c [CM] 



£ := min -, ll , and U := max i -, 1, ~l \ } — 
q } [q [q+l)C 



Since 



< r + {p — qC) X + {1 — C)y for (x, y) e (£, cxd)^ , 
then Cqx + Cy <r + px + qy ioi x,y > C„ i.e., 

r + px + y 



C < 



qx + y 



(x,y) G (0,oo)^ 



From the definition of U we have 



-r + {qU - p)C + {U-l)C>Q 
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Write X, y G [C, oo) as x = C + v, y = C + w ior v, w G [0, oo). Then for v,w G [0, oo), 

-r+{qU -p)x + {l{ -l)y 



= -r + {qU -p){C + v) + {U- 1) (£ + w) 

= -r + {qU-p)C + {U-l)C + {qU-p)v + {U-l)w 

> 

that is, 

< 6/ lor X, 2/ G LL, oo) 

g X + y 

□ 

Inspection of the proof of Proposition [l] reveals that, given that we have Lemma 17 the conclusion of 
the proposition is true concerning Eq.(3-2). The statement is given next. 

Proposition 3 At least one of the following statements is true: 

(A) Every solution to Eq.(3-2) converges to the equilibrium. 

(B) There exist m* , M* with L < m* < M* s.t. 

(i) [m*,M*] is an invariant interval for Eq. (3-2), i.e., /([m*,M*] x [m*,M*]) C [m*,M*]. 

(ii) Every solution to Eq.(3-2) eventually enters [m*,M*]. 

(iii) f{x,y) is coordinate-wise strictly monotonic on [m*,M*]^. 

The proof of Theorem |3] may be reproduced here in its entirety with the only change being the elimination 
of the case r < 0, which presently does not apply. Everything else in the proof applies to Eq.(3-2). The 
proof of Theorem |4] is complete. 
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7 Appendix: Computer Algebra System Code 



(l + x) ( 1 + y) (u^ - u + X + y) 



g;cx_' y_}l = 



(b+l)u^-(g+l)u + X * gy 



y_l = 



b X + y 



T|£x_' y_}l ={f[=<^ yl' =}; 



DELTa2 = g[Cx, y} ] -g[T[T[Cx, y} ] ] ] ; 



z 



ruleb = 



b -» 



z + 1 




/. ruleb] ; 



ruleg = g factor 



s + 1 



g+ 1 



+ "b /. ruleg /. ruleb] ; 



ruleu = u Factor 



b + 1 



nuiiiD2 = Mumer at or [ Together [DELTAS ]] ; 

nuiiilD2vk = Numerator [ Together [ nuiniD2 / . (x-»u + v + k, y-»u + v}]]; 

list lD2vk = Numerator [ Together [ nuinlD2vk / . {ruleu^ ruleg, ruleb}]] /. Plus -> List; 
mln21 = Min[ llstlD2vk / . Cz-»1, s-»l, t-»l,v-»i, k->l}]; 
num2D2wk = Numerator [ Together [ nuiiiiD2 / . {x-»u + v, y-»u+v +k}]]; 

llst2D2wk = Numerator [ Together [ nuin2D2wk / . {ruleu, ruleg, ruleb}]] /. Plus -> List; 

mln22 = Min[ llst2D2wk / . (z-»l, s ^ 1, t-»l,v-»l, k-»l}]; 

num22 = Numerator [ Together [ iiu]iiD2 /. {x->u/(w+ 1), y->u/(v+ 1)}]]; 

num2D2vk = Expand [ nuin2 2 / . {v v + k} ] ; 

ll5t2D2vk = Numerator [ Together [ num2D2vk / . (ruleu, ruleg, ruleb}]] /. Plus -» List; 
mln23 = Min[ llst2D2vk / . Cz-»1, s->l, t-»l,w-»l, k->l}]; 
num2D2wk = Expand [ niun22 / . {w -» v + k} ] ; 

llst2D2wk = Numerator [ Together [ num2D2wk / . (ruleu, ruleg, ruleb}]] /. Plus -> List; 
mln24 = Min[ llst2D2wk / . {z-»l, s->l, t-»l,v->l, k->l}]; 
Print [ "Minijnal Coefficients:", (mln21, mln22, mln23 , mln21}]; 
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[ £x_, y_} 1 = Together [ T[ {x, y} ] /. {b -> , g -» , u -> 1 + s}] ; 

1+s+t 1+s+t+r 

DELTA3 = Together [g[ {X, y} ] - g[ Ts [ Ts [ Tz [ {x, y} ] ] ] ] ] ; 
niunDELTas = Numerator [ step3 ] ; 

nujiiiDEIiTA3 2 = Numerator [ Together [ nuiitDELTA3 /. {x-^l + s + v^ y-^ 1 + s+ w}]]; 
nuiiiDELTa33 = Expand[ nuinDELTA32 ] ; 
nujiiDELTA33vk = Expand [ nuiiiDELTA3 3 /. w v + k] ; 

mln3 1 = Min[ numDELTA33vk / . Plus -» List / . Cs-»l,t->lr r-»l,v-»l, k->l}]; 
nuiiiDELTa33vk = Expand [ nuiiiiDELTA3 3 /. v -» w + k] ; 

min3 2 = Min[ nuniDELTA33wk /. Plus -» List / . {s-»l,t->l, r-»l,w-»l, k-»l}]; 

1 + s 1 + s 

nuiiiDELTA32 = Numerator [ Together [ nuiiiiiDELTA3 /. '[x -> , y -> }] ] ' 

1 + V 1 + w 

nuiiiDELTA3 3 = Expand[ nuinDELTA32 ] ; 
nuiiiDELTA33vk = Expand [ nuiiiiDELTA3 3 /. w -> v + k] ; 

min3 3 = Min[ nuiiiDELTA33vk /. Plus -» List / . {s-»l,t-»l, r-»l,v-»l, k-»l}]; 
nuinDELTA33wk = Expand [ nuiniDELTA3 3 /. v -» w + k] ; 

min3 1 = Min[nuiiiDELTA33vk /. Plus -> List / . Cs->l,t-»l, r->l,w->l, k->l}]; 
Print [ "Minimal Coefficients:", {min31, min32, min33, min31}]; 

Table 2: Mathematica code needed to do the calculations in Claim [3] when g <b The functions g, f and 
T are defined as before (not shown). This input was tested on Mathematica Version 5.0 |19j . 
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